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Abstract 

We investigate primordial non-Gaussianity and dark matter isocurvature fluctu- 
ations in the modulated reheating and the curvaton scenarios. In these scenarios, 
large non-Gaussianity can be generated, on the other hand, depending on how dark 
matter is produced, too large isocurvature fluctuations can also arise, which is in- 
consistent with current observations. In this paper, we study this issue in a mixed 
scenario where the curvature fluctuations can also be produced from the inflaton 
fluctuations as well as those from a light scalar field such as the modulus and the 
curvaton. We show that primordial fluctuations can be highly non-Gaussian with- 
out conflicting the current constraint on isocurvature fluctuations for such mixed 
scenarios. However, if the constraint on isocurvature fluctuations becomes severer 
as expected by the Planck satellite, /nL; a nonlinearity parameter for adiabatic fluc- 
tuations, should be very small as /nl ^ 3, which would give interesting implications 
for the generation mechanism of dark matter. Non-Gaussianity from isocurvature 
fluctuations is also discussed in these scenarios. 



1 Introduction 



Cosmic density fluctuations that we can observe today originate to those generated at 
the early universe. Since current cosmological observations on density fluctuations are so 
precise, they provide us a lot of information on the physics of the early universe. Among 
various observables, non-Gaussianity of primordial perturbations has been attracting much 
attention recently. One of the reason is that current and upcoming cosmological observa- 
tions can well probe the Gaussian nature of perturbations more accurately than before. 
The standard simple inflation model predicts almost Gaussian primordial fluctuations and 
its deviation from the Gaussian fluctuations is less than 10"^. Thus, larger deviation from 
the Gaussian fluctuations, which is still less than 0.1% of the Gaussian part, would indi- 
cate that we need some mechanism of generating primordial curvature perturbations other 
than the standard single slow-roll field inflation. The size of non-Gaussianity is usually 
characterized by a non-linearity parameter /nl- The purely Gaussian fluctuations corre- 
spond to /nl = and the current constraint on this quantity is given as —9 < /nl < 111 
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at 95% C.L. in Ref. [1] and -4 < /nl < 80 at 95% C.L. in Ref. [2 

Another test of primordial fluctuations is the adiabaticity of primordial fluctuations. 
If any deviation from purely adiabatic fluctuations is discovered, it would have important 
implications for the generation mechanism of matter (dark matter (DM) and baryon) in 
the early universe. The adiabatic relation between matter and radiation should be satisfied 
when both of them are created from a single component. If, however, the deviation from 
the adiabatic relation is detected, then it implies that radiation and matter have originated 
from separate components. To quantify the deviation from the adiabaticity, the fraction of 
the isocurvature fluctuations to the total ones is usually used, which is defined with their 
power spectra at some reference scale and denoted as a. Depending on how isocurvature 
fluctuations are generated, such an isocurvature mode can be correlated/uncorrelated with 
the adiabatic ones. The current observational limits on (uncorrelated type) and a_i 
(correlated type) are given by < 0.16 (0.072) at 95% C.L. from WMAP5-only (from 
WMAP5+BA0+SN) and a_i < 0.011 (0.0041) at 95% C.L. from WMAP5-only (from 
WMAP5+BA0+SN), respectively [1]. 

As a possible mechanism of generating large non-Gaussianity, the curvaton [3,4] and the 
modulated reheating scenarios [5] have been investigated in various contexts. Although 
these mechanisms are attractive with regard to producing large non-Gaussianity, when 
one considers the generation of baryon asymmetry and DM in these mechanisms, large 
isocurvature fluctuations may arise [6-11], which would indicate that such scenarios are 
disfavored by cosmological observations. 

In our previous letter [11], we have investigated density fluctuations in a scenario with 
gravitino DM in the framework of modulated reheating [5], which is known to gener- 
ate large non-Gaussianity [12]. Then, we have shown that gravitino DM is disfavored if 
the adiabatic curvature perturbations have large local-type non-Gaussianity because of 



Two types of non-Gaussianity are often discussed in the literatures, one is the so-caUed local type 
and the other is the equilateral type. In this paper, we only consider the local type non-Gaussianity. 
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generating too large DM isocurvature fluctuations simultaneousljCj. However, in general, 
fluctuations from the inflaton can also contribute to cosmic density fluctuations today even 
if we consider mechanisms such as the modulated reheating and the curvaton. In such a 
case, density fluctuations are a mixture of fluctuations originating from multiple sources. 
Such a mixed scenario in the framework of the modulated reheating [14] and the curva- 
ton [15-20] has been discussed by several authors. With the contribution from the inflaton 
fluctuations, it is expected that the situation could be dramatically changed because both 
the isocurvature fluctuations and the non-Gaussianity are suppressed. In this paper, we 
consider the density fluctuations in a scenario with dark matter isocurvature fluctuations 
in the mixed modulated reheating and the mixed curvaton scenarios. Then, we investigate 
whether large non-Gaussianity can be generated without conflicting observational limits 
on isocurvature fluctuations. Although we focus on the gravitino (axino) DM scenario 
when we discuss the issue in the framework of the mixed modulated reheating, our discus- 
sion on the mixed curvaton scenario can apply to a generic DM which originates from the 
inflaton or the curvaton. We also discuss non-Gaussianity from isocurvature fluctuations 
in these scenarios, which can also affect the non-linearity of cosmic microwave background 
(CMB) temperature fluctuations. 

The paper is organized as follows. In section [21 we summarize a formulation to inves- 
tigate non-Gaussianity based on 6N formalism and cold dark matter (CDM) isocurvature 
fluctuations. In section [Sj we briefly review the modulated reheating scenario. In sec- 
tion HI we investigate non-Gaussianity and the gravitino DM isocurvature fluctuations in 
the mixed modulated reheating scenario. In section [5l we move on to the discussion on 
the mixed curvaton scenario. In section [6l we discuss the non-linearity of isocurvature 
fluctuations in these scenarios. Section [7] summarizes our results. Throughout this paper, 
we set the reduced Planck mass M|j = (SttG)"^ to be unity, where G is the gravitational 
constant. 



2 Formulation 

In this section, we summarize a formalism to discuss primordial non-Gaussianity and 
isocurvature fluctuations. 

2.1 Non-linearity parameter in ^A^ formalism 

First, we give the definition of a non-linearity parameter which characterizes non-Gaussianity 
of the primordial curvature perturbations. In order to evaluate the curvature perturba- 
tions ( on super-horizon scales, we adopt the 6N formalism [21]. In this formalism, the 
curvature perturbations on sufficiently large scales ( at the final time t = tj are identical 

'^'^ In fact, in the framework of the curvaton scenario [3,4], which is also known as a good candidate to 
produce large non-Gaussianity [13], such DM scenario may not be viable because of too large isocurvature 
fluctuations as well. 
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to the perturbations of the e-folding number measured in the homogeneous FRW Universe 
from the initial time t = to the final time as 

C(t/)^5iV(t^,t.) , (1) 

where N represents the e-folding number defined as N{tf, t,,) = J^*-^ Hdt with H being the 
Hubble parameter. Usually a final hypersurface at t = is taken to be a uniform energy 
density one and an initial hypersurface a.tt = t^, to be a flat one. Taking the initial time to 
be some time shortly after horizon crossing during inflation, we can expand 6N in terms 
of fluctuations of scalar fields (f"' on the initial flat hypersurface as 

atj) ^ Nj^: + lKb6^:6cpi , (2) 

up to the second order. Here, a superscript a labels a scalar field and Na = dN{tf, t) / dip"'{t) 1^=^^ 
and Nab = d'^N(tf,t)/dip"'(t)dip''(t)\t=t^. The summation is implied for the repeated in- 
dices. 

To discuss non-Gaussianity, one usually considers the bispectrum (3-point correlation 
function) of the curvature perturbations which is written as 

iCkSk^Ck,) = (2^)'^^'^ (^1 + k2 + h) B^ih, h, h). (3) 

To quantify the size of non-Gaussianity, the non-linearity parameter /nl is often adopted 
and is defined as 

B^i^ h, h) = ^/nl (Pc(^i)n(^2) + 2 perms.) . (4) 
In the SN formalism, we can write the 3 point function of the curvature perturbation as 

X {Ps{h)Ps{k2) + 2 perms.) 5^=^) (^i + h + h) , (5) 

where we have neglected the nonlinearity of 6ip'^. Here the indices are lowered and raised 
by using the Kronecker's delta 6°'''. The second term in the parenthesis is a contribution 
from the one-loop correction, km = min{/cj}(z = 1,2,3) and L is a cutoff scale which 
is often taken to be the order of the present Hubble scale. Ps and P,^ represent power 
spectra of fluctuations of scalar fields and the curvature perturbations, respectively. They 
are given and related to Vs and by 

(S^XS^l^^J = {2n)W^ (k, + h) Psiki) = ^S^'Vs6^''> {k + ^2) , (6) 
(C^^C.^^) = (27r)^5(3) {k, + k,) P^{k,) = ^n^^'^ (^1 + ^2) . (7) 



3 



The relation between and Ps can be written as 

P^{k) = [NaN- + K,N-'Vs HkL)] Ps{k) . (8) 
Then we can express the non-hnearity parameter /nl as 

Here, we neglect the contribution coming from the higher order in 5ip^, for example, terms 
with Nabc, N„hrA , and so on. We will justify this assumption for the modulated reheating 



scenario latei*'^ 



2.2 Cold dark matter (CDM) isocurvature fluctuations 

Now let us move on to the issue of CDM isocurvature fluctuations. If we consider the 
curvature perturbations Q on the spatial slices of uniform density pi for the i-th component, 
which are related to the total curvature perturbation ^ as C = ^iCiPi/ Pi isocurvature 
fluctuations between CDM and radiation are defined as 

'S'cDM = 3(CcDM — Cr), (10) 

where Ccdm and Cr are the curvature perturbations defined on the slice of Pcdm and pr 
being uniform, respectively. One can also write down this quantity using the fluctuations 
of the ratio between the number density of CDM, tt-cdm, and the entropy, s, as 

•JCDM — 1 — • U-LJ 

J^CDM/S ncDM S 

To parametrize the contribution from isocurvature fluctuations, one usually uses the 
fraction of isocurvature fluctuations to the total ones which is defined as 

- p^ik,) + Ps{k,y 

where k^ is some reference scale at which the power spectra are evaluated. Ps{k) is the 
power spectrum for isocurvature fluctuations defined by 

{S^,^S^^) = (27r)3<5(3) (^k, + h) Psik,) = ^^5^^'^ (^i + h) . (13) 

In the scenario discussed in the following, isocurvature fluctuations can be correlated/uncorrelated 
with adiabatic ones and both can arise simultaneously. In such a case, we need to define 



For the curvaton scenario, as shown m Ref. [22], one can find that the third order terms with Nabc 
can be neglected for large /nl in the absence of the non-linear evolution of the curvaton field between the 
horizon crossing and the start of curvaton oscillation. 
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the fraction a separately for correlated and uncorrected ones, respectively. Detailed dis- 
cussion on this point will be made in the next section. 

As in the case for the curvature (adiabatic) fluctuations, we can also define the non- 
linearity parameters for isocurvature fluctuations 

'S'cDM = 'S'g + f^^^Sg, (14) 

where Sg is the Gaussian part of isocurvature fluctuations. In some cases, the contribution 
from the second order term dominates over the first order one. In such a case, the frac- 
tion of isocurvature fluctuations a represents the size of non-Gaussianity. Regarding the 
notation for non-Gaussianity from isocurvature fluctuations, we follow those of Ref. [23]. 



3 Modulated reheating scenario 

In this section we give a brief review of the modulated reheating scenario and some re- 
sults for non-Gaussianity of the curvature perturbations and gravitino DM isocurvature 
fluctuations in this scenario. 



3.1 Non-Gaussianity of the curvature perturbations 

Let us start with considering the background dynamics during reheating era, in order to 
evaluate the curvature perturbations in the modulated reheating scenario based on 6N 
formalism. Here we assume that the inflaton oscillates under a quadratic potential after 
inflation and hence the energy density of the inflaton behaves like a matter during its 
oscillation. The homogeneous background equations during reheating era are given by 

r 

+ 3p0 = -— p^, (15) 



dN H' 
^ + 4p. = -p,, (16) 

H' = l (P0 + Pr) , (17) 

where p,^ and pr are energy densities of the inflaton and radiation, respectively. F is the 
decay rate of the inflaton into radiation. In the modulated reheating scenario, the decay 
rate of the inflaton depends on a light scalar field (so-called modulus) a. Thus F can 
fluctuate due to fluctuations of the modulus. The e-folding number N during reheating 
era can be given by 



N{tf,ti)= In 



a(ti) 



where a{t) is the scale factor. The initial time ti and the final time tf are respectively taken 
to be the time at the end of inflation and some time after the completion of reheating. 
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Under the sudden decay approximation, the inflaton decays into the radiation at H = T, 
suddenly. Then, one can rewrite the e-folding number as 

where t = t^ec is the time when H = T. From t = ti to t = tdec, the Universe is 
dominated by the inflaton behaving like matter and then becomes dominated by radiation 
after t = tdcc- Hence, we have 

m^t.) = -h-(T^]-h-(^] , (20) 



3 \H{U)J 2 

where we have used H oc a~^^'^ during matter dominated era and H oc during radiation 
dominated era. In the 6N formalism, the curvature perturbations can be generated from 
fluctuations of the decay rate originating from those fluctuations of the modulus 5cr*, 

C{tf) ~ 5iV = NJa, + ^N^Jal , (21) 



where 



dN 

iV. = —T'ia) , (22) 



dN d^N 9 

N.. = It r'V) + ^ (rV))' • (23) 



From Eq. fl20p . we find 



dN _ 11 

'dT ~ ~6f ' ^ 
d^N _1 1 

9n ~ 6f2 • ^^^^ 

Using these expressions, we can evaluate the power spectrum of the curvature perturba- 
tions and also the non-linearity parameter /nl in the modulated reheating scenario given 
by Eq. m- 

3.2 Gravitino DM isocurvature fluctuations 

Here we discuss isocurvature fluctuations from gravitino DM in the modulated reheat- 



ing scenarig^i. There are two major ways to produce gravitino DM^^I. One is from the 



'^'^ The same discussion also applies to axino DM which are produced by thermal scattering during 
reheating [24]. 

Gravitinos can also be produced by the decay of the next-to-the-lightest supersymmetric particle 
(NLSP). However, the detailed calculations show that the constraints from BBN are severer than those 
from the overclosure of the universe, irrespective of the kind of particle of NLSP as long as it is the MSSM 
particle [25]. Therefore, a scenario with (most) DM being gravitinos can be realized only for the following 
two cases. 
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scattering of particles in the thermal plasma and their relic abundance is evaluated as [26] , 



where Tr is the reheating temperature, ^3/2 is the gravitino mass, m^j is the gaugino 
masses for i-th generation, and gi is the gauge coupling. Notice that the relic abundance 
is proportional to the reheating temperature Tr. The other way is to produce gravitinos 
non-thermally from the decay of some heavy scalar field such as the infiaton or moduli [27]. 
In this case, the yield can be written as 

^ = 1^3/.;^. (27) 

where S3/2 is the branching ratio of the decay into gravitinos. When gravitinos are pro- 
duced from the jets, S3/2 should be understood as those including its multiplicity. Since 
the reheating temperature is related to the decay rate of the infiaton as F oc T^, then 
S3/2 oc T^^, we have the Ti^-dependence of the ^3/2/^, in the case of non-thermal produc- 
tion, as 

^ « ^. (28) 

In the modulated reheating scenario, the reheating temperature fluctuates in space due to 
the fluctuations of the modulus field a, thus gravitino DM isocurvature fluctuations can 
be generated as 

^(nj^^^Wi, 

' "3/2/3 ' ^ ' 

where the positive and negative signs correspond to the cases with thermal and non- 
thermal productions, respectively. Hence, we can easily flnd that if the gravitinos consti- 
tute DM in the Universe, isocurvature fluctuations are generated in both cases. 

From current observations such as CMB, the magnitude of isocurvature fluctuations is 
strongly constrained. Since Tr is proportional to F^/^, we flnd 5Tr/Tr = 5F/(2F). Then 
using the curvature perturbations given by Eq. (!2T]) . 6*3/2 can be related to ( as 

S3/2 = ±^ = ±§ ^ T3C (30) 

Thus we have 83/2/ C — =f3, which is fully correlated to the curvature perturbations and 
already contradicts with the current observations [11]. In fact, too large isocurvature 
fluctuations can also be generated when we consider the production of gravitino DM in 
the curvaton mechanism. Since the modulated reheating and the curvaton scenarios are 
the major mechanisms of generating large non-Gaussianity, if the primordial curvature 
fluctuations are found to be highly non-Gaussian in the future, gravitino DM scenarios 
may be disfavored because of too large isocurvature fluctuations [11]. However, when the 
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curvature fluctuations from tlie inflaton also contribute to today's density fluctuations, 
which we call a mixed scenario, the fraction of isocurvature fluctuations would be diluted. 
Thus in such a case, the DM scenario discussed above may be liberated. On the other hand, 
non-Gaussianity, which can be generated from fluctuations of the modulus or curvaton, 
would also be reduced. Thus it is interesting to investigate how large non-Gaussianity can 
be generated without conflicting the constraint on isocurvature fluctuations in such mixed 
scenarios of the modulated reheating and the curvaton. In the next section, we discuss 
the mixed scenario of the modulated reheating. Then in Section [5l we consider such a 
scenario in the framework of the curvaton. 

4 Mixed modulated reheating scenario 

In this section, we discuss the gravitino DM isocurvature fluctuations and non-Gaussianity 
in the mixed modulated reheating scenario where fluctuations of the inflaton also con- 
tribute to the curvature perturbations as well as those from the modulated reheating. 

4.1 Non-Gaussianity 

In the SN formalism, the curvature perturbations in the mixed case are, up to the second 
order, given by 



where and (^(o-) represent the curvature perturbations originating from the inflaton 
and another scalar field which is assumed to be the modulus here, respectively. They are 
written as 



Here we have neglected the non-linearity coming from the inflaton fluctuations since it is 
of the order of slow-roll parameters and very small. Notice that or Qcr) is different from 
or (cr, which is often seen in the literatures. The former represents the contribution to 
the total curvature perturbation coming from the fluctuations of or o" while the latter 
is the curvature perturbation on the slices of uniform density or po-- Then, the power 
spectrum is given by 



C = C(<^) + C(a) , 



(31) 



(32) 



(33) 



(34) 



where and Pa are power spectra which are defined by 




(35) 
(36) 
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Here we have assumed that and a are uncorrelated. 

For the discussion later, we define the ratio of the power spectra between P^- and 
at some reference scale and express them as 

Nl 



Nl 



■Vs \n{kL) 



with 



Ni 



Nl 



V 

i_ /r 

36 \Y 



1 

Ye 



N„ 



Nl 



6 



r'vr 



(37) 

(38) 
(39) 
(40) 



Here for simplicity we have considered the standard slow-roll inflation model and e is so- 
called a slow-roll parameter. With this definition, the limits of P — > and oo correspond 
to the pure infiaton and pure modulus cases, respectively. By using P, the non-linearity 
parameter /nl given by Eq. ([9]) can be rewritten as 



6. 



P N^ 



[l + RY 



Ni 



(41) 



In principle, the term with N^„„ can also arise in P^^l. However, by a simple inspection 
of such a term, we can see that such higher order derivative terms can be neglected. If we 
take into account the contribution from the third order term in 6a, R can be expressed as 

Nl 
Nl 



R 



Nl 



N„ 



Nl N, 



(42) 



Now we compare Nl^jN^ with NafjfjjN^ and show that the latter is much smaller than 
the former, in particular, when /nl is large. From Eqs. fl39l) and (HOj) . we obtain 



N^ 



.g-3(n(£: 



(43) 



Thus the ratio of these combinations is given by 



N IN 

* (7(7(7 j * cr 

KJNl 



2-3 



rvr 
(Tyif 



1 - 



r'vr 
(rvrf 



-\ 2 



(44) 



In fact, the term with F"' = (ir'Y(a)lda^ can also appear. However, in an exphcit model we consider 
in the following, such higher order derivatives vanish. Thus we neglect such a term. 
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From Eq. (HT]) we can easily find that at least we need N^jfj/N"^ ^ 1 in order to realize 
large /nl and then this condition corresponds to 



r'vr 



> 1 . 



{T'/Vf 

Under this condition Eq. (144|) is approximately given by 



(45) 



N IN 




- Y"/Y - 


GU 1 (7 







< 1 



(46) 



Hence as far as we consider the case with /nl ^ C^(l); the term with N^jcja in the one-loop 
correction is negligible. 

Similarly, we can easily confirm that higher order terms like N„„„„ and N„„„„„, which 
can arise in the one-loop correction term of bispectrum and trispectrum, are also negligible. 



4.2 Gravitino DM isocurvature fluctuations 



Here, we discuss the gravitino DM isocurvature fluctuations in the mixed modulated re- 
heating scenario. Gravitino DM isocurvature fluctuations are given by 



S- 



3/2 



2T' 



-3 



IE! 

2T 



6a! 



1 



5al 



(47) 



which are correlated with the adiabatic fluctuations as S'3/2 — — 3Co- as discussed in the 
previous sectional. Since here we consider a mixed scenario where the curvature fluc- 
tuations can also be generated from the inflaton, only some fraction of the isocurvature 
perturbations is correlated with the (total) curvature perturbations. Then, S'3/2 can be 
divided into the following two parts. 



S'3/2 — S'corr + S'l 



Sr. 



corr 

Prr 



3R 



"3-^ (Cw + Cm) = -yTr^' 

SR A P^. \ 3 

TTr l^^^^ " tJ^"^ 



R 



{RCw - Cm) 



(4^ 



"^"^ In fact, the relation between the isocurvature and adiabatic fluctuations slightly deviates from S'3/2 = 
— 3Ccr due to the non-linear terms. However, in order to realize large non-Gaussianity, the condition 
given by Eq. ([IS]) must be satisfied. Hence, by neglecting N^/Ncra in the non-linear term of S'3/2 we 
approximately obtain S3/2 — — 3Ccr- 
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where the first part is denoted as S'corr and the second as S'uncorr- Using the expression ([3 
we have the following relations, 

R 1 



1 + R ' 1 + R 



To express the size of the contribution from isocurvature fiuctuations, we define the ratio 
of the power spectrum relative to the total one as 

ftcorr = p , p , (50) 



■^corr 



Ps 

-June 



ncorr / CT 1 \ 

"uncorr = p , p ' \^^) 

uncorr 

for the correlated and uncorrelated parts, respectively. Following the notation of Ref . [1] , 
we define the cross-correlation coefficient /? as 

(3 = , (52) 

where Psco^^c denotes the cross-correlation power spectrum defined by 

(^corr.^.C,^.) = (2vr)'5(') (^1 + h) P5_c(^i) ■ (53) 

Hence the correlation coefficient is (3 = -|-1(— 1) for the thermally (non-thermally) 
produced gravitino DM in the modulated reheating scenario. 

As already mentioned, the size of isocurvature fiuctuations is now severely constrained 
by observations of CMB and so on. In fact, in the present model, correlated and uncorre- 
lated isocurvature perturbations coexist, thus we need to take into account both contribu- 
tions simultaneously to obtain observational constraints on Ocorr and auncorr- However, such 
analysis is not available in the literatures. Hence as reference values, we adopt the con- 
straints on Ocorr and auncorr obtained separately from recent WMAP5 results: Ocorr < 0.011 
and auncorr < 0.16 at 95 % C.L. from the WMAP-only analysis [1@. 



4.3 A Simple Model 

Now let us work on some explicit model of the modulated reheating scenario. We consider 
the following interaction between an infiaton </> and a fermion ip: 

C,nt ~ -9{<y)H^ , (54) 

In fact, the constraint for correlated isocurvature fluctuations here is obtained for (3 = —1. However, 
it is expected that the sign of the correlation does not affect the constraint on the size much. Thus we 
refer this value regardless of the sign of /3 as a reference value. 
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where the couphng constant g depends on a modulus field a. We assume that g{a) can 
be written as 



9o 



(55) 



where M is some energy scale, gi and g2 are some coefficients and \<t\/M < 1 is assumed. 
The decay rate of the inffaton through this interaction is given by 



which implies that the form of the decay rate of inffaton is 



a 



M 



(56) 



(57) 



where A and B are some coefficients. Substituting Eq. (1571) to Eqs. (139|) and (HOjl . we 
obtain 



Ni 



1 



36M2 



A + 2B (cj/M) 



and 



6M2 



Vl + A{a/M) + B{a/My 
A + 2B (a/M) 



{51 



1 + A{a/M) + B{(T/My 



2B 



1 + A {a/M) + B {a /My 



(59) 



Now we discuss the non-linearity parameter /nl and the size of isocurvature fiuctua- 
tions in this model. There are four parameters in the model; M, cr/M, A and B. Here, 
we will show the results for /nl and the isocurvature fraction by considering the limit of 
\a\/M < 1 or A ^ 0. 

In the case of |(j| /M <C 1, in Eq. and N^^ in Eq. (ISUil are respectively reduced 



as 



Nl 



A 

6M 



JSl ~ 



{A^ - 2B) 



On the other hand, for the limit of A 

-\ 2 



6M2 

0, these are written as 



(60) 



Nl 



1 



36M2 



2fi((7/M) 
+ B {a/Mf 



N ~ 



6M2 



25 (a/M) 
1 + 5 (a/M)^ 



25 



1 + B {a /My 

(61) 

In Fig. [1], we show contours of /nl along with ttcorr (left panel) and ttuncon- (right panel) 
in the M-A plane. We have fixed the value oi B a,s B = —1.0 and assumed |cr| /M <C 1. 
In Fig. [21 we show the same but in the M-|i?| plane. In the figure, we have fixed as 
A = 1.0 and assumed \a\ /M <C 1. In Fig. El we have fixed the value of model parameters 
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Figure 1: (color online) In the left (right) panel, contours of /nl and Ocorr (^uncorr) in the 
M-A plane are shown. We have fixed as i? = —1.0 and assumed |cr| /M <^ 1. Red thin 
lines in the left panel show contours of Ocorr = 0.011 (solid line) and 0.0041 (dashed line). 
In the right panel, the lines are shown for Ouncorr = 0.16 (solid line), 0.072 (dashed line) 
and 0.01 (dotted line). Green shaded regions are constrained by the current observational 
limit for the isocurvature fractions from WMAP5-only analysis. Blue thick lines show 
contours of /nl = 5 (dotted line), 10 (dashed line) and 100 (solid line). 




Figure 2: (color online) In the left (right) panel, contours of /nl and Ocorr (ttuncorr) in the 
M-\B\ plane are shown. We have fixed as A = 1.0 and assumed |(t| /M -C 1. 
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Figure 3: (color online) In the left (right) panel, contours of /nl and Ocorr (ttuncorr) in the 
M-|(t| /M plane are shown. We have fixed as B = —1.0 and A = 0. 

as B = —1.0 and A = 0, then show contours of /nl as well as Ocorr and auncon- in the 
M-|(t| /M plane. Regarding the slow-roll parameter e for the infiaton, we take e = 0.01 
in the following analyses. In all figures, red thin lines in the left panel correspond to 
contours of Ocorr = 0.011 (solid line) and 0.0041 (dashed line). These numbers correspond 
to the 95 % C.L. limit from WMAP5 only and WMAP5+BA0+SN analyses, respectively. 
In the right panel, they correspond to contours of auncorr = 0.16 (solid line) and 0.072 
(dashed line). These numbers again correspond to the 95 % C.L. limit from WMAP5 only 
and WMAP5+BA0+SN analyses, respectively. In the near future, we will have a more 
stringent limit from PLANCK [28-30] and its projected limit on uncorrelated isocurvature 
fluctuations will be Ouncorr < 0.01, thus we also show the contour of ttuncorr = 0.01 with 
dotted line. Green shaded regions are constrained by the current observational limit for 
the isocurvature fractions from WMAP-only analysis. Blue thick lines show contours of 
/nl = 5 (dotted line), 10 (dashed line) and 100 (solid line). 

From these figures, we can see that large values of /nl ~ (10 — 100) can still be realized 
without conflicting the current constraint on isocurvature fluctuations. In particular, when 
the one-loop correction term dominates, the following simple relation holds between /nl 

and Ouncorr- 

p o / '^uncorr \ "^^^ /■nr,\ 

In Figs. [1] and [3l the regions where the one-loop term dominates correspond to downside 
of the figures and in Fig. [21 it corresponds to upside of the figure. In these parameter 
regions, the present limit on the isocurvature fluctuations still allows relatively large non- 
Gaussianity in the scenario. 

Here it should be mentioned that the isocurvature constraint is severer for uncorrelated 
one. The reason is as follows. Here we consider the case where the adiabatic curvature per- 
turbations are mainly generated from fluctuations of the infiaton field, which corresponds 
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to the case with R <^ 1. On the other hand, the gravitino DM isocurvature fluctuations 
are produced from fluctuations of the modulus field as given in Eq. (1471) . Hence, in the 
small R region, the gravitino DM isocurvature fluctuations have mostly become uncor- 
related type. This fact can be also understood by noting Eq. (H9l) . In the near future, 
the data from PLANCK will be available and the constraint on isocurvature fluctuations 
would be much severer than that of the current one. The projected limit on uncorrelated 
isocurvature fluctuations would be a^ncow < 0.01 at 95 % C.L. [28]. Thus, for reference, 
we also plot the corresponding contour in the figures, from which we can find that, if 
we obtain such a stringent constraint, large non-Gaussianity cannot be generated even in 
the mixed scenario. Thus if in the future, large local-type non-Gaussianity is confirmed, 
but the isocurvature constraint becomes as severe as that mentioned above, gravitino DM 
scenario would be ruled out even if we consider the mixed scenario. 



4.4 Brief comments on the trispectrum 

Before closing this section, we would like to make brief comments on the trispectrum in 
the mixed modulated reheating scenario. The local type trispectrum can be parameterized 
with two non-linearity parameters (7nl and tnl as [31] 



Tnl {Pdki3)Pdh)Pdh) + 11 perms.) 
54 

+ ^9NL iPdk2)Pdks)Pdh) + 3 perms.) 



5(3) (^1 + k2 + h + h) 



(63) 



where k. 



ki -\- kj 



scenario we have 



By adopting the 5N formalism, in the mixed modulated reheating 



R 



[l + Rf 



(64) 



including the one-loop correction. From Eqs. fHTl) and ([6; 

l + R 



Tnl 



R 



fi 



2 

NL • 



(65) 



Hence, for small R there is a possibility of generating large tnl which may be detected 
in the future experiments, without contradicting the current observational constraint on 
the isocurvature fluctuations. As discussed above, once we obtain a severe constraint 
on isocurvature fluctuations, the value of /nl would be small in the mixed modulated 
reheating scenario. However, it is possible that a non-Gaussian signature comes from the 
trispectrum but not from the bispectrum. 
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5 Mixed curvaton scenario 



Now in this section, we consider CDM isocurvature fluctuations in the curvaton scenario. 
Isocurvature fluctuations in the curvaton have been investigated in the literatures [6- 
10]. Here we investigate this issue in the framework of a mixed scenario and focusing 
on how large non-Gaussianity can be produced without conflicting with the constraints 
on isocurvature fluctuations. Notice that the discussions given in this section apply not 
only to gravitino (axino) DM but also a generic CDM, although we focused on such DM 
candidates in the previous section. As pointed out in [8,11], when the number density of 
CDM freezes before the curvaton decay, too large isocurvature fluctuations are generated 
and excluded by cosmological observations. However, this conclusion is valid for the 
original curvaton scenario in which fluctuations from the curvaton are only responsible for 
density fluctuations today. But, in general, fluctuations of the inflaton also contribute in 
addition to those from the curvaton. Such a mixed scenario has been extensively studied 
for the adiabatic fluctuations in [16-20] and for baryon isocurvature fluctuations [10]. Here 
we consider CDM isocurvature fluctuations in the mixed curvaton scenario for the cases 
with CDM being produced from the decay of the inflaton and/or the curvaton, paying 
particular attention to how large non-Gaussianity can be without conflicting with the 
isocurvature constraint. 

In the same way as in the mixed modulated scenario, under the sudden decay ap- 
proximation, the adiabatic curvature perturbations on the uniform (total) energy density 
hypersurface are analytically given by [20] 

C = C(<^) + C(<7) 

C(.) = NJa. + N^Jal = ^/dec^ + ^/dcc (l - ^/dec - lfl?j , (66) 

where and a denote the inflaton and the curvaton, respectively, and a is taken to be 



positive without loss of generalitjl^. We have neglected the non-linear part of fluctuations 



of the inflaton since they are very small. Here /dec is defined by 

3p<x 



/ 



dec 



4pr + 3Pa 



(67) 



where and pr are respectively energy densities of the curvaton field and radiation. 
t = tdec is the time at the curvaton decay. For large non-Gaussianity, at least we need 



3 1/ 4 2 2 

' 1 - ^/dcc - -/dec > 1 , (6^ 



4/decV 3^"^^ 3' 

which leads to /dec ^ 1- Hereinafter, we consider the case where this condition is satisfied. 



The case where the second term dominates the hnear term in Q^.-) corresponds to the "ungaussiton" 
scenario [32-34]. Our analysis includes such a scenario automatically. 
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After the curvaton decay, the CDM isocurvature fluctuations are expressed as 



'S'cDM — 3 (CcDM — C) ) 



(69) 



where C is given by Eq. (166|) and Ccdm is the curvature perturbation on the uniform CDM 
energy density hypersurface. 

Regarding the production of CDM, one can consider two cases: dominant residual 
CDM is generated from the decay of the inflaton or the curvaton. In the following, we 
consider each case separately. 

5.1 Case with CDM from the inflaton decay 

First, let us consider the case where dominant residual CDM is generated from the decay 
of inflaton, which implies that CDM has the same fluctuations of the inflaton, that is, 
Ccdm = C<f> which denotes the curvature perturbation on uniform inflaton energy density 
hypersurface. In the case where /dec <^ 1, we can consider ^ C(9i) hence Ccdm — C(0)- 
Then, we obtain [20] 



Since we are considering a mixed scenario, the isocurvature fluctuations here should have 
uncorrelated and correlated parts with adiabatic fluctuations which originate from both 
the inflaton and the curvaton. With the deflnition of the correlation coefficient of Eq. fl52l) . 
P = 1 for this case. Then we separate isocurvature fluctuations into two parts as done in 
the previous section. 

In Fig. m we plot contours of /nl as well as Ocorr (left panel) and Ouncorr (right panel) 
in the cr,,-/dcc plane. The results are quite similar to those for the mixed modulated 
reheating scenario discussed in the previous section. At the current level of the constraint 
on isocurvature fluctuations, relatively large non-Gaussianity can be generated without 
conflicting the constraint. However, once the limit becomes severe as ttuncorr < 0.01 which 
is expected in PLANCK, /nl should not be large in this case. 

5.2 Case with CDM from the curvaton decay 

Next, let us consider the case where dominant CDM component is generated from the decay 
of the curvaton, which implies that Ccdm = Cai^" ({a} = /dccCo-) because /dec <^ 1, where Co- 
is the curvature perturbation on the uniform curvaton energy density hypersurface. Then, 
we can obtain [20] 




(70) 




(71) 
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Figure 4: (color online) Contours of /nl along with Ocorr (left panel) and ttuncorr (right 
panel) in the cr^-fdcc plane. Red thin lines in the left panel show contours of Ocorr = 0.011 
(solid line) and 0.0041 (dashed line). In the right panel, we show contours of auncorr = 
0.16 (solid line), 0.072 (dashed line) and 0.01 (dotted line). Green shaded regions are 
constrained by the current observational limit on isocurvature fractions from WMAP5 
only. Blue thick lines are contours of /nl = 5 (dotted line), 10 (dashed line) and 100 (solid 
line). 



In the case where /dec 1, Sqdm is approximately given by 



CDM 



with 



/dec 
2 /dec 



(72) 



2/< 



dec 



0-* 3 

from Eq. (IMj) . Then, the power spectrum of the isocurvature fluctuations is given by 

1 



(73) 



PscUk) ^ 9jT- + NLVs HkL)] Ps{k) 

J dec 
J dec 



(74) 



In fact, in this case, we can easily find that CDM isocurvature fluctuations become large 
even for /nl = 0(1). From the above equation, the ratio of the power spectrum of the 
CDM isocurvature fluctuations to that of the adiabatic (curvature) perturbations is given 
by 



P' 



S'CDM 



9 



_9_Pa ^ 

•I dec C -I dec 



R . 



(75) 
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From Eqs. (137|) . we obtain the following inequality: 



Nl 9' 



R>^ = T^e^ , (76) 



where we have used = (2e) ^ and Eq. fl73l) . Hence, we have the inequality given by 



> 8^ . (77) 



On the other hand, from Eqs. (14T!) . assuming that i? -C 1, the non-linearity parameter 
/nl can be written as, 

N 4 e 

/nl ^ = -UecR- . (78) 

3 at 

Combining Eq. (l78l) with Eq. ( 1771) . we obtain 

> R-'f^,lhL. (79) 



RScOM 



When <C 1 and /dec < 1, the fraction of the power spectrum of the CDM isocurvature 
fluctuations to that of the adiabatic (curvature) perturbations becomes much larger than 
the size of the non-Gaussianity, Pscom/Rc ^ /nl- Thus, in the case where dominant 
CDM component is generated from the decay of curvaton, large non-Gaussianity from the 
(adiabatic) curvature fluctuations cannot be generated without conflicting observational 
constraint on CDM isocurvature fluctuations even at the current level. 



6 Non-Gaussianity from DM isocurvature fluctuations 

So far we have discussed non-Gaussianity of adiabatic fluctuations. However, isocurvature 
fluctuations can also produce large non-Gaussianity, which has been discussed recently 
in [20,23,35,36]. Thus in this section, we investigate non-Gaussianity from isocurvature 
fluctuations in our scenario. 

First, we consider non-Gaussianity from isocurvature fluctuations in the modulated 
reheating scenario. As discussed in Sec. HI the isocurvature fluctuations can be written, 
up to the second order, as 



3/2 



ir 
2T 



6a ^ 



1£! 
2T 



5at 



(80) 



Depending on the explicit form of F and the value of (5(T*, the discussion on non-Gaussianity 
can be divided into two cases. When the linear term dominates over the second order one in 
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Eq. (IHOl) . it is convenient to define a non-linearity parameter for isocurvature fluctuations 

/■(iso 
'NL 



as in the case of adiabatic ones. We denote the non-linearity parameter as f^^^ and define 



it as the same as the counterpart in the adiabatic case: 

'S'3/2 = Sg + f^^^Sg, (81) 

where Sg is the Gaussian part of isocurvature fluctuations. With this definition, is 
calculated as 

^(iso) _ r /r ^ 1 

/nl - ^^,^^^2 2 ■ ' 

Notice that f^^^ takes almost the same value as /nl for adiabatic fluctuation although the 
sign is different. However, even if the size of is comparable to that of the counterpart 
for the adiabatic mode, the signal from the bispectrum depends on the combination of 
a^/]^L°^ and af^^^ for the bispectrum coming from the 3-point function of (SSS) and 
{SS(), i.e., purely isocurvature and correlated parts, respectively. Thus they are sup- 
pressed by the fraction of isocurvature perturbations a, which is severely constrained by 
observations. Hence non-Gaussianity from isocurvature fluctuations cannot be large com- 
pared to the adiabatic one in this case. 

However, when the second order term dominates over the first order one in Eq. ( IHOll . 
the argument becomes different. In this kind of case, the isocurvature fluctuations can 
be simply written as 5*3/2 = Sg, which is called the quadratic model in [23]. (The model 
characterized by Eq. (IHTj) is called linear model in [23].) In the quadratic model, its power 
spectrum is determined by the second order term, thus non-Gaussianity in this model 
can be represented only with 03/2, which characterizes the size of the power spectrum of 
isocurvature fluctuations relative to the adiabatic ones and is given by 

«3/2 - -p^^- (83) 

This fraction should be small to fit to current data of CMB power spectrum, thus ^ 
PS3/2 ^^'^ using Eq. ( H9l) and i? ^ 1, we obtain the following relation, 

(1 + R)P, + Ps,,, -l + R P, - ^'^^ 

As discussed in the previous section, this parameter is constrained as auncorr < 0.072 at 
95 % C.L. Thus we use this value as a representative one in the following. 

Here it should be noted that, even if we characterize the primordial non-Gaussianity 
for both adiabatic and isocurvature fluctuations with the nonlinearity parameters /nl and 
f'^^\ their evolutions to the present epoch, which are encoded in the transfer functions, 
are different. Thus the comparison of the non-linearity between isocurvature and adiabatic 
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fluctuations is not so simple. However, some useful relations are obtained among /nl, /nl 
and a-i/2 in [23]. Thus we make use of those relations. 

Since non-Gaussianity from isocurvature fluctuations tends to be small for the case with 
the linear model, as discussed above, we here consider the quadratic case. Non-Gaussianity 
from the uncorrelated term can be estimated effectively as 



where /nl in the left hand side is the non-linearity parameter for adiabatic fluctuations. 
This relation has been obtained by finding 03/2 which gives the same S/N for adiabatic /nl 
assuming WMAP5 noise. As seen from the above relation, relatively large non-Gaussian 
fluctuations can be generated from isocurvature fluctuations even if we take the value 
allowed by current severe constraint on the isocurvature fluctuations. 

Next we discuss non-Gaussianity from isocurvature fluctuations in the curvaton sce- 
nario. For the curvaton scenario, we considered two cases for the generation of CDM and 
discussed isocurvature fluctuations for each separetely. When CDM is created from the 
inflaton, isocurvature fluctuations are given by Eq. (170|) . from which we can evaluate the 
non-linearity parameter /nl°'' 



Notice that f^i^ is almost the same in size as the adiabatic counterpart /nl except from 
the sign, which is similar to the case of the modulated reheating scenario. Although the 
size of the non-linearity parameters are almost the same between /nl and /nl""*, the signal 
of the bispectrum is suppressed by the isocurvature fraction as discussed above. Thus 
non-Gaussianity from the isocurvature fluctuations in this case would be also small as 
well. 

The other case we considered is that CDM is produced from the decay of the curvaton. 
In this case, the isocurvature fluctuations are written as Eq. (171]) . then the nonlinearity 
parameter can be given by 



As seen from this expression, when /dec is close to 1, /nl° would be very large. Thus in this 
case, even if the bispectrum itself is suppressed by the fraction of isocurvature fluctuations, 
its signal can be very large. However, notice that the sign of /nl''* negative. In fact, when 
/dec is close to 1, /nl for the adiabatic fluctuations is /nl ~ C'(l). Thus, when /dec ~ 1, 
non-linearity mainly comes from isocurvature fluctuations and its size can be large but the 
sign is negative. Since too large value of /nl°'' would be disfavored by observations, this 
scenario may contradict even with current observations. 
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When the quadratic term dominated over the first order term, the argument is the 
same as the case with modulated reheating, in which the fraction of isocurvature fluc- 
tuation a gives the size of non-Gaussianity or bispectrum. Hence, when the quadratic 
term dominates in C(o-) , non-Gaussianity from isocurvature fluctuations can be large in the 
curvaton scenario as well. 



7 Summary and Discussion 

In this paper, we considered DM isocurvature fluctuations and non-Gaussianity in models 
where the adiabatic curvature fluctuations can be produced not only from of a light scalar 
field other than inflaton (modulus or curvaton), but from the inflaton fluctuations, which 
is called a mixed scenario. Regarding the non-Gaussianity of the curvature (adiabatic) 
fluctuations, we have found that relatively large non-Gaussianity can be realized as /nl ~ 
(9(10 — 100) without conflicting with the current constraints on the fraction of the CDM 
isocurvature fluctuations. In other words, the current limit on isocurvature fluctuations is 
not severe enough to prohibit large /nl for such mixed scenarios. However, for the future 
CMB experiments such as Planck satellite, the limit on the uncorrelated isocurvature 
fluctuations will be improved as auncorr ^ 0.01 [28-30]. We showed that this projected 
limit translates into the bound on the non-linearity parameter as /nl < 3. 

In fact, although we have mainly discussed non-Gaussianity in the curvature (adiabatic) 
fluctuations, nonlinearity can also arise from isocurvature fluctuations. We have also 
discussed non-Gaussianity of this type and found that it can be large as shown in Eq. fl85l) . 
which corresponds to the size of the adiabatic nonlinearity parameter as /nl ~ 10 even 
with the projected Planck limit for Ouncorr- However, it should be noted that the signature 
in the bispectrum of adiabatic and isocurvature fluctuations are not th e sa me, thus we 



may differentiate non-Gaussianity from these fluctuations to some extenl i^^i . Thus, if we 
find that the value of the local type (adiabatic) non-linearity parameter is large in the 
future experiments, cosmological scenarios with gravitino DM may be disfavored even if 
we consider a mixed fluctuation scenario. 

Furthermore, we have also investigated DM isocurvature fluctuations in the framework 
of the curvaton, in particular focusing on a mixed scenario. Our discussions for the curva- 
ton case also apply to a generic DM, which originates from the inflaton or the curvaton. 
We have found that large non-Gaussianity in the curvature (adiabatic) perturbations is 
possible with the current level of isocurvature constraints in the curvaton well. 
However, as in the case of the modulated reheating scenario, if the limit becomes severer 
and large (adiabatic) non-Gaussianity of the local type is found in the future, DM are 
unlikely to be produced from the decay of the inflaton or the curvaton in the curvaton 
scenario, which would give important implications to the generation mechanism of DM. 



"^^^ The phase difference between adiabatic and isocurvature fluctuations in acoustic oscillations can help 
to distinguish non-Gaussinity from these two modes [37]. 
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